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Abstract 



We describe flux tubes and their interactions in a low energy sigma model 
induced by SU(Nj) — > SO(Nf) flavor symmetry breaking in SO(N c ) QCD. 
Unlike standard QCD, this model allows gauge confinement to manifest itself 
in the low energy theory, which has unscreened spinor color sources and global 
Z2 flux tubes. We construct the flux tubes and show how they mediate the 
confinement of spinor sources. We further examine the flux tubes' quantum 
stability, spectrum and interactions. We find that flux tubes are Alice strings, 
despite ambiguities in defining parallel transport. Furthermore, twisted loops 
of flux tube support skyrmion number, just as gauged Alice strings form 
loops that support monopole charge. This model, while phenomenologically 
nonviable, thus affords a perspective on both the dynamics of confinement 
and on subtleties which arise for global Alice strings. 
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I. INTRODUCTION AND MODEL 



Nonlinear sigma models have had great success in describing low energy QCD phe- 
nomenology. However, they have not captured the hallmark feature of QCD: confine- 
ment, where the potential between qq pairs grows linearly with separation. Such lin- 
ear potentials can arise as tension-carrying flux tubes - - between q and q, for Yang- 
Mills QCD, or between any unscreened sources which persist in low energy effective the- 
ories. The conventional Skyrme model, induced by the global flavor symmetry breaking 
SU(Nf) L x SU(Nf) R — > SU(Nf) diag , supports neither unscreened sources nor flux tubes. 
Global flux tubes are classified by 7i2(G/H), due to their constant vev at spatial infinity. 
This necessarily vanishes whenever the vacuum manifold G/H is itself a Lie group, as occurs 
in the Skyrme model. 

Witten noted that an SO(N c ) gauge theory of QCD induces a different Skyrme model, 
whose topology can support flux tubes. This theory has Nf lefthanded quarks qi, which 
transform as (real) fundamentals under color SO(N c ). SO(N c ) has no gauge anomaly, and 
thus requires no right-handed quarks; however, its U(1)b anomaly breaks baryon number 
to Z 2 . The high energy theory thus has the full symmetry group Z 2 x SO(N c ) x SU(Nf). 
At low energies, flavor SU {Nf) breaks to its SO(Nf) subgroup, due to formation of a quark 
condensate (qf^ Ey qlj) (where a and i are color and flavor indices respectively). This 
condensate interacts with quark excitations in the theory, inducing an effective Majorana 
mass for the quarks 



which displays the SU(Nj) — > SO(Nf) flavor symmetry breaking explicitly. The Gold- 
stone modes E are described by a SU(Nf)/ SO(Nf) nonlinear sigma model, with skyrmions 
(tt 3 (G/H) = Z 2 , for Nf > 4) and flux tubes (ji 2 (G / H) = Z 2 , for Nf > 3). This defect 
classification makes physical sense: for baryons, identified with antibaryons because real 
quarks are identified with antiquarks; and for flux tubes, whose Z 2 structure emerges as a 




(1.1) 



2 



response to external spinor sources, which can be screened by fundamental quarks only in 
even combinations. 

In this paper, we elaborate on our recent work with Manohar 0, constructing the flux 
tubes in this theory and showing that their interactions with skyrmions and spinor sources 
obey heuristic expectations. Confinement of spinor sources in an SO(N c ) gauge theory can 
thus manifest itself in the low energy sigma model, through relic phenomenology in the 
presence of unscreened sources. 

We organize our results as follows: section II derives the unique flux tube form with 
minimal energy; section III examines its classical stability and dynamics; section IV its 
quantum stability and spectrum; and section V its interactions and relationships with other 
fundamental objects in the theory. The early sections reveal that minimal flux tubes lie in a 
planar subspace of the vacuum manifold SU(Nf) / SO(Nf) ; while section V discusses why our 
flux tubes are Alice strings, despite ambiguities in defining parallel transport around them; 
why they carry skyrmion number when twisted; and how they can be viewed as mediating 
the confinement of spinor sources. 

II. FINDING NONTRIVIAL FLUX TUBES 

To construct topological defects, we seek configurations where the condensate £ varies 
spatially in both nonsingular and nontrivial ways. Nonsingular variations assume a form 
dictated by the transformation properties of S. Demonstrating their nontriviality, however, 
is complicated in our modified Skyrme model. The difficulty stems from the same source as 
the flux tubes themselves: the fact that the global vacuum manifold G/H is not itself a Lie 
group, but only a quotient space. Thus, unlike skyrmions in the standard Skyrme model, 
winding numbers for E cannot be obtained from a group structure on G/H alone, as indices 
dependent only on S. Instead they must be determined by homotopy arguments from the 
embedding of £ in G. 

The SU{2) L x SU{2) R Skyrme model presents a simpler case. Here symmetry group 
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transformations qr) act on G by left multiplication and leave £ invariant under the 
diagonal subgroup. £ can thus be constructed from underlying group transformations: £ = 
<?lW j transforming as £ — > LluR) under the group element (L,R). When E = 1 , L — R 
gives the diagonal unbroken subgroup. For S^l, however, the embedding H C G changes 
and L = R = g rotates E on the vacuum manifold, E — > <?£<?L Using these transformation 
laws, and imposing spherical symmetry to minimize energy, we may write an arbitrary 7r 3 
defect as E = L(r2)E G (r)i? t (f2). We choose our basis to diagonalize S (r) = exp {iF(r)r z }, 
with F(r — > oo) = 0. E then has finite energy only when L = R, keeping E constant at 
spatial infinity, and is singular at the origin unless F(0) = nn. 

Establishing this nonsingular form E = g(Q) exp {iF(r)r z }g^(Q) requires only the 
quotient space structure G/H = (SU(2) L x SU(2) R ) /SU(2) diag . We require the full Lie 
group structure of G/H, however, to show that a particular g(Q) induces nontrivial 7r 3 , 
or equivalently, to reduce the Pontryagin index to a function of E. Since ir^G/H) = 
n 3 (SU(2)) = Z, the homotopy index 7r 3 simply measures the number of times E(r, Q) cov- 
ers SU(2). This is easily determined by group integration or inspection. For the skyrmion, 
g(Q) = exp {— i0T z /2} exp {— i<fyr y /2} rotates r z to an arbitrary Lie algebra element f • r. 
Thus E = exp {iF(r)r ■ f} consists of the exponentiation of the full Lie algebra from to 
nit — by definition, covering the Lie group n times, for a 7r 3 winding n. 

Our model yields a nonsingular ansatz for E as directly as the standard one. G = 
SU(Nf), acting on itself by left multiplication, leaves E invariant under the orthogonal 
subgroup. We can thus write E as gg T , transforming as E — > aEa T under the group element 
a. When E = 11, a G SO(Nf) gives the unbroken subgroup H; elsewhere, the embedding 
H C G is parallel transported to gHg^ . Again, a cylindrically symmetric defect can be 
written as an r-dependent vev, with angle- dependent group rotation: E = g(0)Yi o (r)g T (9). 
Choosing E c = 11 at spatial infinity restricts g(6) to lie in the unbroken SO(n) subgroup, 
and to commute with E D (r = 0), if E is to be nonsingular with finite energy. 

We can further characterize E in terms of su(Nf) basis generators. These specify ro- 
tations in all 2-dimensional subplanes (jk), and are usually taken as follows (with Cartan 
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norm tr T a T b = \8 ab ). The rank Nf — 1 Cartan subalgebra has as basis diagonal matrices 
T d , 

T d = (2d(d + 1) Y 1 ' 2 diag(l, . . . , 1, -d, 0, . . . , 0), (2.1) 

with ones in the first d entries and d — 1 to Nf — 1. We can overspecify this basis, sacri- 
ficing the Cartan norm, by the set {\T z (jk)} in an subplanes of the Nf dimensional vector 
space. Off-diagonal generators {^r x Qk) , \ T y(jk) } complete the su(Nf) basis. Flavor sym- 
metry breaking divides the basis into two sets: antisymmetric matrices {Th} = {\r y (jk) }, 
generating H = SO(Nf); and symmetric matrices {T b } = {\r x (jk)i ^ T z(jk)}, generating no 
symmetries. S (r), a unitary symmetric matrix, can be written exp {iF b (r)Tb}. Thus S as- 
sumes the most general nonsingular form E = exp {iF b (r)T b )}h~ 1 (9), for /i(6*) e and 
F b {r) ranging from zero at infinity to 2nn 5 bb r (for some fixed direction b') at the origin. 

Having obtained nonsingular configurations E, we must demonstrate their nontriviality. 
Unlike the conventional Skyrme model, this requires an understanding of how E = gg T arises 
from the underlying group mapping g. For flux tubes, we must construct E from the exact 
sequence 

tt 2 (SU(Nf)) = -> tt 2 (SU(Nf)/SO(N f )) 

-> n 1 (SO(N f )) = Z 2 -> tti (SU(N f )) = . 

That is, gg T gives a nontrivial E only if g corresponds to some mapping from the plane to 
SU(Nf), with boundary values in the SO(Nf) subgroup. Furthermore, when parametrized 
as a family of loops, g must start at the identity and end on a nontrivial loop in SO(Nf). 
Taking (a G [0, 27r], /5 G [0, 7r]) as our coordinates on the plane, these criteria become 

!U when a = 0, a = 2tt, oy 8 = 
(2.2) 
h 2 (a) when /3 = 7r, 

where h 2 (a) is a nontrivial loop in SO(Nf). Such nontrivial loops can be written h 2 (a) = 
exp {ia 2nhT h }, where T h is the set {^T y ^k)} introduced above and UhT h generates rotations 
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in a single plane. (Of course, these loops can be deformed, but deformations from geodesic 
form lengthen the loop and induce additional gradient energy in E. We discard them, to 
focus on mappings that produce minimal gradient energy.) 

We now show that the embeddings g(a, (3) have minimal gradient energy only when they 
induce flux tubes S of a unique form. We show this by imposing consistency and minimal 
energy conditions on the most general trivialization g(a,P). We then relate the topological 
coordinates (a, (3) to cylindrical coordinates (r, 0), to obtain the physical flux tube S(r, 9). 
Finally we consider the low-lying deformations of S(r, 6) which can be favored by potential 
energy terms. We thus obtain a family of non-trivial flux tubes, among whom a minimal 
representative is selected dynamically. 

We construct the most general trivialization g(a,(3) as follows. Start, at P — 0, by left 
multiplying h(a) by its inverse: g{(3 = 0) = h~ 1 (a) h(a) — 1. As (3 varies, allow the left 
multiplier h~ 1 (a) to vary over the full group G: g(a, (3) = gi(P) h^ 1 (a) g 2 (P) h(a). The 
a = condition and uniqueness of the inverse imply ^ 2 (/?) = 9\ 1 {P)- To probe structure on 
G/H (where a — > aa T identifies cosets {gH}), write g\{(3) = b((3) h((3), for h{(3) G H and b 
generated by broken generators T b . (Specifically, since g\{(3) gj((3) = exp {iG b (/3)T b }, we can 
write g\{(3) = exp {iG b (/3)T b /2)} h(f3).) Thus we may write the most general trivialization 
of the loop h 2 (a) in G: 

g (a,(3) = b{(3) h{(3) h-\a) h-\P) 6 _1 (/3) h(a). (2.3) 

This induces minimal energy only in geodesic form 

b{(3) = exp {iip n b T b } 
h((3) = exp {imp h h T h } 

h(a) = exp {ia n h T h }, (2.4) 

where n b and rih,hh are unit vectors over the range of T b ,Th respectively. 

We now show that consistency reduces the distinct choices of b, h, and h. For simplicity, 
take h — 1 at first. Choose a basis in su(Nf) so that n b T b = |r 2 ( 12 ). Then, under 
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conjugation by b(/3), 2nhTh = n^k) T y (jk) rotates as follows: 

b(f3) 2n h T h fo" 1 ^) =n(i2) (cos//3r„(i2) + sin If} t x ( 12 )) + ^ r y(i*0 

2<j<fc 

+ H n 0'fc) (cos(//3/2) Ty^-jfej - (-1) J sin(Z/3/2) r x(ifc) ) . (2.5) 

3=1.2 
fc>2 

To obtain g(a,{3 — ir) — h 2 (a), this conjugated generator must give —2nhT h at (3 — ir. This 
can occur for only two lowest winding possibilities: I can be 1, with tiuja = en 2) vanishing 
outside the (12) plane; or I can be 2, with ny*.) describing a plane that intersects (12) in 
a single line. The second possibility is further constrained by the boundary condition at 
a = 2ir, which requires h(2n) to commute with n^Tb. This occurs only when [n^Tfe, (n^T^) 2 ] 
vanishes — that is, when the plane nu^) intersects (12) along a coordinate axis. Thus only 
three distinct candidates arise for the pair b,h: I = 1, n^-) = e(i 2 ); / = 2,n(y) = e(i 3 ); and 
I = 2,n,(ij) = e( 23 ). 

Allowing nontrivial h(/3) produces no additional flux tubes of minimal energy, as we 



demonstrate explicitly in an appendix. Thus the form [2.3| for the flux tubes's embedding 
reduces to two candidates: 

g (a,p) = b{p) hr l {a) b-\l3) h(a), (2.6) 

with 

!exp {i{3 T bn } 
(2-7) 
exp {i2f3 T bn >} 

where □ and □' denote planes intersecting along a coordinate axis, and T b ,T h are |r 2 , ~r y 
in the indicated planes. Written in this form, we see that the two candidates are in fact 



the same. From equation |2.5| , both choices for b{(3) induce the same rotated group element 
b(/3) h~ 1 {a) & _1 (/3) (modulo sign redefinition of (3 ). We can thus choose a single deformation 
g(a,/3), exploring only a planar SU(2) subgroup of SU(Nf), to produce the flux tube of 
minimal gradient energy associated with the Z2 loop h 2 (a). Note that g indeed gives a Z2 
object, with g(a,/3 = 2ti) = 1. (This Z 2 structure appears in the (3 — > —f3 equivalences 
above, as well). 



From this form we construct the flux tube E = gg T : 

E = b(J3) hr\a) b- 2 {(3) h{a) b(J3). (2.8) 

This gives 

S(a,/3) = 11 - 2 sin 2 (a/2) sin 2 /3 11 □ + zsin 2 (a/2) sin (2/3) r zD -isinasin/? r xD , (2.9) 

where 11 □ gives the identity in the plane □ and vanishes outside it, and 11 is the usual 
SU(Nf) identity. Of course, this form for E can deform while remaining a nontrivial flux 
tube. In particular, for E on G/H, the variables a and (3 give coordinates on the physical 
plane R 2 , identified to S 2 by the condition E — » 11 at spatial infinity (i. e. on the boundary 
a G {0,27r}; /3 G {0,7r}). We can deform E to a radially symmetric form by identifying 
r = oo with this boundary, and r = with the center a = n, (3 = it /2 of the deformation g 
producing E. This yields 

E(r,0) = /i(0) 6(r) ^ _1 (0), for 

(2.10) 

/i(0) = exp{i6» T^n}, 6(r) = exp {«F(r) T bn }, 

with T^n, Tft D as in our definition |2/?| for ^(a, (3). 

We derive this radially symmetric form for E and fix F(r = 0) as follows. We expand 



equation 2.10 to obtain 



E(r, 9) = 1 + ( cos(F/2) - 1 ) 11 □ + i sin (F/2) ( cos 9 r zD - sin t x0 ) . (2.11) 



which can be identified with equation p.9|, term by term. This gives an undeformed relation 
between the topological description g(a,/3) and the spatial form F(r),9: 



F(r) = 4 tan" 1 ( [sin(a/2) sin/3]- 2 - 1 ) 

= 2tan" 1 (cot (a/2) cos/3 + sin/3( [sin(a/2) sin/3]" 2 - 1 V ' ~ 



(2.12) 



where tan 1 gives values in a single period [— 7r/2, vt/2]. The arguments to tan 1 range from 
to ±oo for 9, and from to oo for F(r). Thus 9 covers a range [— 7r,7r], suitable for a 
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polar angle, and F(r) ranges from to 2tt. F(r) assumes the extremes of its range at the 
boundary r = oo, where F(r) = 0, and the center r = 0, where F(r) = 2n. These boundary 
conditions comprise the only relevant feature of the mapping, the homotopy invariant F(r = 
0) — F(r = oo). They establish boundary conditions for any single-winding flux tube of form 
F(0) = 2vr, F(r -> oo) = 0. 



Finally, we note that while form |2.10| for £ indeed minimizes gradient energy, poten- 
tial terms for £ — like a quark mass term — can favor a vev other than exp (iF(r) T^n) ■ 
Potential terms consistent with unbroken SO(Nf) symmetry give a vev we can always di- 
agonalize. We thus consider the extending £ so that F(r) T& n — > Fd(r) Td, where F\ = F 
and Td varies over the Cartan subalgebra generators [2.1| . The nonplanar T^>i commute with 
both T bn and T ha ; that is, they commute with the full embedding g(a,/3). We may thus 
obtain our extension by the simple modification g(a,/3) — > g(a,j3) exp {iF d> ±(r) T d>1 /2}. 
Taking i^>i(r) to vanish as r — > oo, this does not affect the behavior of g on the boundary 
- hence it leaves the topology unchanged. It changes the resulting flux tube only by the 
desired overall multiplication S — > S exp {ii^ >1 (r) Td > ±}. Thus the true minimum energy 
flux tube is one of a family of nontrivial configurations: 

E(r,0) = h{6) b{r) h-\6), with 
h(6) = exp{t6 T hD }, b(r) = exp {tF d (r)T d }, (2.13) 
and Fi(r = 0) = 2tt; F d {r -»• oo) = 0. 
Which of these candidates is realized remains a question of dynamics. 

III. FLUX TUBE STABILITY AND DYNAMICS 

Studying flux tube dynamics begins with the question of stability. As minimal model for 
the Goldstone field E, with stable skyrmions, we have the Skyrme lagrangian 

Co = S tr d^d^ + J-, tr ^d v Y] 2 . (3.1) 



9 



Unaltered, this lagrangian implies a size instability for all flux tubes. For under the rescaling 
E(r, 9) — > E(Ar, the tension of a finite flux tube can always decrease. Specifically, its 
quadratic contribution stays invariant, while its quartic rescales by a factor A 2 — leading to 
an energy minimized when A = 0. Physically, this corresponds to flux tubes which diffuse to 
infinite size to lower their energy. 

To stabilize the flux tubes, we must consider modifications of the minimal forms for E 
and C Q . First, we note that adding higher derivative gradient terms to C Q , analogous to the 
Skyrme term, cannot both stabilize E and produce a positive definite Hamiltonian. This 
holds because such terms can be at most second order in time derivatives, and, by Lorentz 
symmetry, r and 9 derivatives; thus they can never give tension contributions scaling more 
strongly than A 2 . This leaves only two options for stabilizing E. First, we can consider a 
potential for E, giving a tension component that scales as A -2 . Second, we can allow the 
flux tube to vary along its axis — giving E some z-dependence — to obtain a quadratic 
contribution to the energy scaling as A -2 . The second option, implemented by exciting zero 
modes along the flux tube's z-axis, has a structure paralleling that which arises in quantizing 
the flux tube. We thus defer a study of its dynamics until the next section on quantization, 
and focus on stabilization by a potential. Ultimately both approaches give similar results: 
both give as minimal flux tubes representatives from the family [2.13| with vanishing Fd>i(r). 
Stabilization by z-rotation has a richer dynamical structure, however — partly because it 
introduces unconstrained rotational parameters into the problem. 

A stabilizing potential, on the other hand, arises naturally in our theory. By giving 
the quarks bare Majorana masses, we explicitly break the SU(Nf) flavor symmetry of our 
original QCD gauge theory. This induces a pion mass term, 



in the limit of degenerate quark masses. The mass term stabilizes the flux tube while 
affecting its dynamics in a simple way, as we discuss below. 

To study dynamics, we calculate the energy density of flux tubes of the general form 




(3.2) 



10 



2.13|. The Skyrme action EO gives a gradient contribution 



16 



-j-ti{(F' d T d y + - (T 2 n -[F' d T dl T D ] 



(3.3) 



where 



f n =b-\r) [T hQ , 6(r)] . 



(3.4) 



and r has been rescaled into the dimensionless units eF^Tphys- This form for p follows 
algebraically, using no information about the generators T d and Thn- T n measures the 
noncommutativity of the radial generator F d (r) T d with its angular counterpart Thn- The 
nonplanar T d> \ commute with all generators Thn, T d — making T n and [ F' d T d , T a ] indepen- 
dent of -Frf>i. This has two consequences: first, nonplanar terms contribute only the positive 
definite sum J2d>i F d 2 /2 to the gradient energy density. Second, the T n -dependent terms 
span only a planar SU{2) subgroup, where T$=\ = \t z u and Thn = \ T yn- They are thus 
calculable, giving the flux tube gradient energy density 



Po 
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| {|&^ 2 + ^ (1 - cos^) (l + ^' 2 )} 



(3.5) 



The pion mass term |3.2| also contributes to the energy. Form |2.13j for the flux tube yields 
its contribution in closed form: 

,2 



Pr, 



ml 



16e 2 



2N f - 4 cos £ u d/2 cos(Fi/2) - 2 ^ cosuo d 



(3.6) 



xd>l 



d>l 



with 



-d F, 



UJ d 



E 



(3.7) 



>2d(d+l) i >d ^2i(i + l) 
where d — 2, . . . , Nf — 1 and r has again been rescaled to eF^ r p h ys - 

The minima of this potential fix boundary conditions for £ at spatial infinity. Such 
minima have F\ = 2nm, consonant with our a priori boundary condition F\(r — > oo) = 0. 
To minimize with respect to the nonplanar u d , we set 
dp„ 



2 sin [ ^d/2 ] cos(Fi/2) + 2 sinu d = 

Kd>l 



(3.8) 
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with positive second derivatives d 2 p m /du!^, d 2 p m /dF 2 . This gives the degenerate family ujd = 
2irjd, for any integers jd summing to the same parity as m. As the simplest representative, 
we take u>2 = 2nm and Ud = for d > 2. This completes our boundary conditions at infinite 
radius: F\ — > 0, Ud — > 0, for all d; corresponding to Fd —>■ for all d, in our original variables. 

Away from spatial infinity, we must propagate these boundary conditions inward to 
obtain flux tubes obeying the full equations of motion. However, given our boundary con- 
ditions ujd(r — > oo) = 0, these equations trivialize for the nonplanar variables Fd>i- They 



have independent gradient energy terms |3.5| , which combine with the potential qJ^ to give 
equations of motion 

£ 1 B, {rl® = ford>l. (3.9) 

16 r oF d oui e 



From equation |3.8|, the vacuum values uj e = 2irj e cause the source terms to vanish on the 



right hand side — even when the planar field F± departs from its vacuum value. Thus our 
boundary conditions, setting Fd — uJd = at infinite radius, induce only constant solutions 



as we propagate equation |3.9| inward to the origin. This gives Fd that vanish identically for 
nonplanar d > 1. 

This implies that the minimum energy flux tube for any Nf varies only over the planar 



SU(2) subgroup. It takes the simple form |2.10| , with F(r) = -Fi(r), and has energy density 

P = § [\F' 2 + i (1 - cosF) (l + F' 2 ) + A 2 (1 - cos(F/2))} , (3.10) 

where A = 2m n /eF 7T . This determines a nonlinear equation of motion for F: 

( l + i s in 2 (F/2)) - d r (rF') = 
V r z J r 

sin(F/2)| ^cos(F/2) ( 1 - F' 2 ) + ^ s in(F/2) + . (3.11) 



Numerical solutions to equation |3.11| are shown in Figure 1 for different values of A, 
including the physical Ao = 0.236 (e = 2tt, m n = 138 Mev and F n = 186 Mev). Increasing A 
raises the flux tube's energy density while shrinking its core size. The asymptotic regimes 
agree with limits to equation |3.11| . Inside the core, F falls linearly from its boundary value 
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of 27r at the origin; outside, it scales as the hyperbolic Bessel function F ~ x~ l l 2 exp {— x}, 
for x = Ar/2. The parameter 2 A -1 thus sets the flux tube's core size in dimensionless units. 
In physical units, this gives core size m" 1 , sensible for a bound state of Goldstone modes S. 
It carries tension proportional to F 2 — which is proportional to the number of colors iV c in 
a large N c QCD limit. This supports its interpretation as a mediator of the confining force 
between spinor sources, whose scales as the spinor Casimir N c . Numerically we find tension 
4.6.FJ when A = Aq. As A varies, the tension varies as shown in figure 2. 



IV. QUANTIZING THE FLUX TUBE 

A. Semiclassical Zero Modes and Quantum Stability 

Under quantization, the flux tube samples not only the ground state above, but also 
its zero modes. Recall, from section II, E's transformation laws: S — > aSa T under global 
rotations a, giving the residual symmetry H = SO(Nf) at infinite radius where £ = 1. 
However, because S varies nontrivially, global H rotations are not symmetries within the 
core. Instead they produce distinct degenerate configurations, coincident at spatial infinity. 
The zero modes of £ explore these configurations: 

Z(t,r,9)=A(t)Z(r,9)A- 1 (t), (4.1) 

where A(t) = exp {ieF n u t n^T/j} rotates in H with dimensionless frequency uj. f] These 
zero modes have rotational energy confined to the string core; an energy that can be made 
arbitrarily small, classically, by taking uj —>■ 0. 

Classically, we calculate this rotational energy from the Skyrme action O: 

P. = 4^tr (f 2 -[F' d T d ,f\ 2 --\T D ,fA . (4.2) 



16 V L a 1 J r 



*We neglect other excitations, such as bending modes. 
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Here T n , from equation [3.4| , measures the noncommutativity of radial and angular generators 
F d {r)T d and T ha , while 



f = b-\r) \ h-\6) n h T h h(8) , b(r) } 



(4.3) 



measures the noncommutativity of A(t) with E(r, 9). We have retreated to the general form 



2.13| — including nonplanar contributions — for E(r, 9); thus we retain complete expressions 



375| and |3l] for its gradient and potential energy. We do this for two reasons: first, we must 

insure that the planar vacuum F d> \ = — favored by the weak pion mass term — remains 

minimal despite quantum fluctuations due to zero modes; and second, we wish to explore 

stabilization by ^-rotation, where A = A(z). 

Calculating the rotational energy |4T2] is straightforward. We expand 2n/ l T/ l as n^ k j T y (j k y, 

then note that the tilde operation, defined in equation |3.4| , acts on basis elements as follows: 
/ 

—t{iijk — 1 ) 

'cos uj jk - 1) T y(jk) + smuj jk T x(jk) . (4.4) 



Here iTj k is given by 



n jk = exp {iujj k } = exp {iF d (T dt k -T djj )} , 



(4.5) 



where T d> k denotes the fc-th diagonal entry in the Cartan generator T d , from equation |2TT 
The commutator terms act on the relevant basis elements r x (jk), T y (jk) to give 



[ T d , r a (j k ) ] = -i{T dj k - T dj j) e ab T b (j k ) Vj, k 
[ 2T n , r a n> ] = ip (cosu; n - 1) r a n» + i sinu; n e ab r b u" 



(4.6) 



where e ab is the two dimensional permutation matrix e xy . The second commutator contributes 
only when planes □ and □' intersect, with [T y n,T yn >] = ip T yD » defining □" and p = ±1. 
Altogether, these expressions imply the rotational energy 



T?2. ,2 

7T 
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TV 



(12) 



(1 - cosiq) (l + F{ 2 ) + £ nf jk) (1 - coso; ifc ) (l + u' jk 2 ) + 



k>j>2 



k>2 



2^ n {lk)+ n 



2 ^ 
(2k), 



' (1 -cosuj lk )(l + uj' lk 2 + ^(1- COS Fi)) ^ 



i +(1 - coscusfc) (l+uj' 2k 2 + ^2 (1 - cosFi 
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/ J 



(4.7) 



This form for the rotational energy depends on the planar Fx, both explicitly and through 
the variables Uj k . To extract its Fx dependence, we note that 



where 



Ujk — — Wj-l, 



(4.8) 



id>l / \rf>i 



(4.9) 



and the nonplanar u>d>i are given by eq. pT7| . We may thus rewrite the rotational energy as 



FW 



Pu 



n 



(12) 



;i-cosF 1 )(l + F 1 ' 2 )+ £ 



n 



k>j>2 



COS 
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n 



fc>2 



(lk) 



+ n 



(2k) 1 



V 



/ J 



(4.10) 



■l + S' k _ 1 2 + F{ 2 /4 + ^(l-cosF 1 ) 
• (1 — cos Sk-i cos(Fi/2)) 
+F{ S 1 ^ sin S k _x sin(F 1 /2) 

where S k = Wfc + | Sd>i ^rf- This depends only on Fi (explicitly) and the nonplanar c<Jd>i 
(explicitly and through Sk-i, Wjk)- 

Neglecting the pion mass term for now, our boundary conditions require F\ and to 
approach a minimum of the potential [4.10 at spatial infinity. As in the previous section, 
these minima occur when Fx = lum and u d>1 = 27rj d , for any integers jd whose sum has 
the same parity as m. We take as the simplest representative uj% = Im-n and uJd>2 = 0. (The 
specific case Nf = 3 has additional degenerate vacua, non-coincident with those of the pion 
mass term, given by uj 2 = 2m7r/3; however, these reduce to the pion mass vacua above when 
more families are added.) Thus our boundary condition Fx(r — > 00) = induces the full 
boundary conditions at infinite radius: Fx — > 0, u>d>x —* 0; corresponding to F d — > for all 
d, in our original variables. 

The rotational energy generates spatially-dependent corrections to the equations of 
motion for Fd- However, these corrections trivialize for the nonplanar variables i*d>i, given 



our boundary conditions uj d = 0. We show this as follows: the rotational energy |4.10 
introduces an additional term 
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P- (^--dr(r^)) (4.11) 
dF d [du K r r \ duj'J K } 



into the right hand side of the equation of motion |3]9| for F&. This term is generally quite 
complicated; however, when the nonplanar u>d>i assume their asymptotic values of zero, it 
simplifies. For a single uj e (before multiplying by the Jacobian), it becomes 



u 2 



E e+ i + \Y, E k \ , with 



k>2 



E k = -F[ S'^ 8^(^/2) - i d r ( 2rS' k _ 1 ) (1 - cos(F 1 /2)) . (4.12) 



Asymptotically, where Fi — 0, all A e vanish, and the full equation of motion |3]9] is again 
sourceless. This gives asymptotic solutions for the nonplanar Fd which are again Bessel 
functions of order zero — implying, for finite energy, that Fd>i and all its derivatives (or 
u>d>i and all its derivatives) vanish at infinite radius. 

However, given nonplanar terms that obey u>d>i = oj' d> \ = ^>i = 0> the perturbations 
A e vanish, regardless of the planar field F\(r). So, as for pion mass stabilization, equation 
[3]9| remains sourceless as we propagate it inward to the origin, giving F d> ± that vanish 
identically. 

Thus the planar vacuum F d> ± = survives quantum fluctuations due to zero modes. 



Its static limit retains the form [2.10| , with F(r) = i*i(r) varying only over a planar SU(2) 



subgroup. It has classical energy, from equation [4.10| , 



Fl \F' 2 + (r- 2 + u 2 nl) (1 - cosF) (l + F' 2 ) + A 2 (1 - cos(F/2)) 
Ptot - — S } (4.1d) 

+u 2 n 2 D , (1 - cos(F/2)) (l + F' 2 /4 + ^ (1 - cos F)) 
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where 



nh = n {X2) , n n , = J2 k>2 n (lk) + n (2k) , (4.14) 

describe the orientation of the zero mode rotation A(t) relative to S(r, 6). 

For a flux tube stabilized by z-rotation, where A = A(z), the rotational terms in p to t 
deform the nonlinear equation of motion for F . Classically, we can rescale uj to set n 2 +n 2 , = 
1 (a rescaling unnecessary when Nf = 3). This gives energy as a monotonic function of n 2 - 
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which measures the component of £ 's z-dependence due to slip dislocation, where E(z, r, 9) 
rotates about its core in a single internal space plane □, but starts at the x-axis with 
changing offset angle 9q(z). Its converse, n^,, measures the twist dislocation in Yl(z,r,$): 
the extent to which A(z) rotates the internal space plane, □ — > R a b (z)D, in which £(2, r, 9) 
cycles about its core. 

Numerical solutions for the z-stabilized F are shown in Figures 3 and 4. Their asymptotic 
behavior agrees with analytic limits: inside the core, F falls linearly from its boundary 
value of 2tt at the origin; while outside, it scales as the hyperbolic Bessel function F ~ 
x -1 / 2 exp{— x}, for x = (A 2 + uj 2 (1 + Sn 2 )) 1 ^ 2 r/2. That is, the flux tube becomes more 
compact as either uj or n a grows. As in the mass stabilization case, shrinking core sizes 
correlate with growing tensions. Thus we see, from Figure 3, that the value n n = is 
favored, and twisting dislocations cost less energy than their slip-offset counterparts. Figure 
4 confirms the flux tube's tendency both to shrink and to gain energy with increasing 
rotational frequency, and shows how significant the rotational deformation of F(r) is. Thus 
the tension not only acquires rotational energy terms, 

t -> t + \uo 2 (n 2 D A n +n 2 D , A n ,)) ; (4.15) 

it also contains hidden rotational dependence through deformations of the ground state 
tension r and moments of inertia A n ,A n /. We explore these deformations in Figure 5, 
showing that r grows linearly with u, due to flux tube compression. The moments of inertia 
instead fall rapidly, dropping by a factor of five as uj grows from to 2, before stabilizing 
at roughly A n = 10 F 2 ,A D ' = 4F 2 for uj > 2 . We show below that this range uj > 2 is 
typical of quantum rotational excitations of the flux tube; thus we are justified in neglecting 
rotational deformation of A n , A D / over this range. 

Finally, we note that the geodesic parametrization assumed above for the rotation A(t), 
while useful for discussing the flux tube's classical limit, does not restrict our analysis. Any 
function A(t) over H = SO(Nf) induces the energy 4.13 , with ujuu^a defined by 

un m = -iieF^y 1 tr A ] A r y{jk) . (4.16) 
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B. The Quantized Spectrum 



From the previous section, there are two ways to stabilize the flux tube: by exploiting 
the pion mass term in the Lagrangian or by z-rotation. These give the two classical solutions 
E(r, 9) (equation |2.13|) and E(r, 9, z) (equation |4.1| , with t — > z). Both have rotational zero 



modes induced by A(t) e H. We now quantize the spectrum of these zero modes to find the 
flux tube's quantum numbers and lowlying excitations |J|. 

Quantizing E(r, 9, z) is a complicated task, because it involves two rotations, the z and t 
rotations. Furthermore, there is no physical input for the z-rotation frequency uj. Therefore, 
we restrict ourselves to quantizing the mass-stabilized solution E(r, 9). 

Equations |4.15| and |4.16| determine a two-dimensional Lagrangian for the flux tube: 



L{z,t) = -f - ^(tr^iT^) 2 - (tr AT m ,) 2 , (4.17) 

where f is a ground state tension and z = eF^z measures dimensionless length along the 



flux tube. The moments of inertia A n and A D / come from equation [4.13 
A n = [ rdr (1 + F' 2 )(l - cosF) 

^ d ' = ^fI ^ ( 1 + F ' 2 /4+ 2^(1 -cos F)) (l-cos(F/2)) . (4.18) 
We note that form [4. 11\ for the Lagrangian comes from only two facts: from the confine- 
ment of E to the planar SU(2) n subgroup in internal space; and from cylindrical symmetry 
in physical space. The skyrmion in this theory, which we discuss in the following section, 
shares these characteristics. Thus its one-dimensional Lagrangian L(t) also has form 4.17 , 
with f, A n and A n / dependent upon skyrmion dynamics. 

To quantize such a Lagrangian, we must do two things: first, we must understand the 
transformation properties of E under the symmetries of L; and second, we must write L in 
terms of invariants of the quantized Noether charges of those symmetries. The first task is 
facilitated by considering the most general form of E on the vacuum manifold: 

Efat) = A(t) e iF{x * )nb{Xl)Tb A-\t) 

— e iF{xi)n b {xi)R b , b T b i ^ (4.19) 



where Rb'b{t) is an orthogonal matrix encoding the isospin rotation due to /1-conjugation: 

A ] TiA = R tj Tj Rijit) = 2 tr ( A\t) T, t A(t) Tj ) . (4.20) 

While Rij is defined for all SU(N) generators T i; it breaks into block diagonal form, 
Rhh' over the space of unbroken and broken generators respectively. 

For the general form E(xj, t), the Lagrangian 4. 17| involves a general sum — 4p(tr A^ATh) 2 . 



This Lagrangian is invariant under global h G SO(Nf), implemented by either left or right 
multiplications A — > hA or A — > Ah. For infinitesimal transformations h e = 1 — 2iehTh, 
S(xj, t) transforms as: 

5 L A = {-i2e h T h )A S L E{xi,t) = -ie h [ 2T h , E( Xi ,t) } 

S R A = A{-i2e h T h ) 6 R E{xi,t) = -ie h A [ 2T h , E(^) ] A^ . (4.21) 

Therefore the left transformation corresponds to SO(Nf) isospin rotation, and its Noether 
charge satisfies an so(Nf) algebra after quantization. The right transformation corre- 
sponds to SO(Nf) isospin rotation in the body fixed frame, with Noether charge T' h . To 
understand the physical interpretation of I' h , we calculate the commutator in the second line 
of equation |4.21| , obtaining 



A[2T h , X(xi)]Ai = 2if hW n b d nv E(x i ,t) . (4.22) 

For the specific cases of flux tube and skyrmion, n b (xi) has two special properties: n b = 
n bn , lying only in an SU(2) subplane of SU(Nf); and n b (xi) is spatially axisymmetric, 
depending linearly on spatial direction components u x and u y . These properties identify 
planar isospin rotations with spatial z-rotations. We see this in two steps: the planarity of 
n b reduces equations |4.21| and |4.22j , when h = hO, to 



5 R E(r,n b (r,Ui)) = -2e f hDb n b >n n ba d nyn E(r,n' b (r,Ui)) , (4.23) 

which the Jacobian relating n b to Ui, at fixed u z , equates with spatial ^-rotation : 

5 R E(r,tt») = -2e e zij u { d Uj £(r,i/») . (4.24) 
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Thus X' D , which generates planar body-centered isospin rotations R, must equal the defect's 
spin operator J z , generating spatial ^-rotations. Other isospin rotations, which destroy the 
planarity of n&, are distinct from other spatial rotations, which change the defect's axis of 
axisymmetry. 

Under right and left infinitesimal transformations of A(t), the rotations Rf/b, which give 
physical coordinates for £(xj,£) on the vacuum manifold, transform as 

SljRb'b = —^hfhb'cRcb 

5nRb'b = —ZthRb'cfhcb ■ (4-25) 

That is, Rb'b is left or right multiplied by {h e ) bc = 1 + ieh(2if hbc ). This shows how isospin 
rotation acting directly on the space of broken generators can obey an so(Nf) algebra - 
specifically the so(Nf) algebra which occurs as a subalgebra of su(Nf) in the adjoint repre- 
sentation, restricted to the block- diagonal component 2ifhb c - 

To find Zfr and X' h , we compute the variation of the Lagrangian under the two infinitesimal 
transformations: 

5 L L = e h X h and 5 R L = e h l' h . (4.26) 

This gives 

X' h = i~K h tx{A^AT h ) (no sum), and X h = Rw^, . (4.27) 

We can now write the defect Hamiltonian in terms of the physical quantized Noether 
charges: the isospin X and angular momentum J . Combining the Lagrangian |4.17| with 
Noether charges ^4.27| yields a Hamiltonian dependent on X' 2 and X^,, where □' are planes 
intersecting □ in only one line. However, the defect S(r, Ui) is invariant under a subgroup 
SO(Nf — 2) of SO(Nf), acting on the subspace Dj_ orthogonal to □. This implies two things: 
X' 2 = I' 2 + J2n>Z'n>] and the defect's allowed quantum states X' must contain a singlet under 
the SO(Nf — 2) subgroup. For Nf > 3, this excludes spinor representations — constraining 
the skyrmions and flux tubes to be bosonic. Lastly, because X is simply a rotation of X', 
and X' D = J Z i we write the Hamiltonian as 
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H{z,t) = r + -ij 2 2 + -J-( X 2 - j; 2 ) , (4.28) 
2A 2A / 

subject to the constraints above: X' n = J Z1 with bosonic defects for Nf > 3. 

Using the numerical results of Section IVA, we may be more explicit about the flux 
tube's quantized energy levels. Relating the rotational energies of equation [4.28| to their 
counterparts in the tension [4. X5| (for dimensionful z) gives 



u 2 = (eig 4 ( j; 2 /4A D 2 + ( X 2 - J 2 )/4A n , 2 ) « 4j; 2 + 25( X 2 - j; 2 ) . (4.29) 

Thus u assumes values to > 2 for quantized zero modes. Over this range, An and An' vary 
little but f is described by the linear fit f = (9.0 + 3.7 (u — 2)) F n /e, as shown in Figure 5. 
Thus our model predicts 



H(z, t) = ^ (l.Q + 3.7 \/4J z 2 + 25( X 2 - J 2 ) + 20 J z 2 + 50( X 2 - J 2 ) (4.30) 
for the allowed excited states (X, J z \ and 

if(S,f)=4.6^ (4.31) 
for the ground state of the flux tube. 



V. FLUX TUBE INTERACTIONS 

We now consider how flux tubes interact with fundamental objects in our low energy 
theory. First, we establish that flux tubes can be Alice strings — forcing some charges and 
particle wave functions to become double-valued. Like gauge theories with Alice strings 
and monopoles, our sigma model has twisted flux loops which form point defects. We con- 
struct these defects — the skyrmions — and demonstrate their nontrivial topology. We then 
consider how they interact with the flux tubes, showing that only topologically trivial com- 
binations of two flux tubes can end on skyrmions. This suggests the physical interpretation 
that, while the baryons in this theory are not confined, the spinor sources which combine to 
form them are, with confinement mediated by the Z2 flux tubes joining them. 
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First we consider twisted flux tubes. We note that our flux tubes may share the defining 
property of Alice strings, which arise in the symmetry-breaking of certain gauge theories. Q 
Alice strings have unbroken symmetries which preserve a local vev but cannot be extended 
globally, since they become multivalued when parallel transported around the spatially vary- 
ing vev of the string. This "Alice" nature is not topologically invariant, but depends on the 
specific Wilson line integral producing the string's asymptotic winding. Whether our flux 
tubes have this trait is similarly ambiguous: using the physical embedding |2.10| of the flux 



tube, we can take g{r,6) = h(9) b(F(r)/2). This gives parallel transported unbroken gen- 
erators g Thw that are double-valued in 9 for each radius. However, all generators can 
be made single-valued by choosing instead g(r,8) = h{9) b(F(r)/2) h(9), which induces the 
same flux tube E. 

Because £ can be viewed as an Alice string, we might expect it to share a property of 
gauged Alice strings in models with monopoles: that twisted string loops support monopole 
charge. For gauged Alice strings, this is necessary, as monopoles alter their charge when 
passing through string loops; moreover, specific examples indicate that twisted Alice loops 
can carry a single unit of monopole charge. ]5|] For global Alice strings, the suggestion arises 
by analogy. In our case, it gains further support from the hidden gauge character of the 
low energy global theory. Specifically, we can recast the flavor-dependent quark mass £ in 
the usual way, as an interaction between flavor gauge fields g _1 (r, 9) g(r,9) and shifted 
quarks Ql = g(r,9) q^- This accomplishes two goals: it equates skyrmion number with 
anomalous baryon current, while identifying the flux tube with a gauged string, whose 
Wilson line integral is Pexp^A-dl However, the Alice ambiguity persists: choosing 
g(r,9) double-valued gives a gauge theory with Wilson loop U(2n) = 1 — 2 H n , making 
both the quarks Ql and the generators T hn / double- valued. A single- valued choice for g(r, 9) 
produces instead single-valued quarks and generators, with Wilson loop U (2n) = 1 . 

This ambiguity can be resolved physically, by considering adiabatic transport of quarks 
around the flux tube. |]] Under such transport, the quarks remain in their mass eigenstates. 
These are governed by two terms: a flavor-independent bare Majorana mass M (breaking 
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SU(Nf) — > SO(Nf) explicitly) and the flavor- dependent Majorana mass /iS, where /i is 
the SU(Nf) -breaking vev. The two determine mass eigenstates that are double-valued in 
9. However, the states have mass splitting Am 2 = — 2iM(/i — /i*) sin(F/2), where we have 
chosen M real. Thus the quarks are degenerate, and unaffected by transport around the 
string, unless the vev \i is misaligned in phase with the bare mass M. In that case, quarks 
at finite radius have double-valued wave functions and Aharonov-Bohm scatter off the flux 
tube, which we identify as an Alice string. [ 

Motivated by its gauge theoretic counterpart, we calculate the skyrmion number of a 
twisted flux tube: 

E(z, r, 9) = A(z) S(r, 9) A~\z) , (5.1) 

where A(z) = exp (iz I n/jT/j) 6 H. Imposing 27r-periodicity in z constrains I: for planar 
rih = no, I must be integral; for nonplanar rih, it must be even. Thus any twisted flux 
loop can deform to the planar flux loop H(z,r, 9) = h(9 + Iz) b(F(r)) h^ 1 (9 + Iz), which 
takes values only in the two-flavor subspace SU(2)/ SO(2) ~ S 2 . This simplification allows 
us to identify the loop's 7r 3 index with its Hopf number, i.e. the linking number between 
any two fibers of constant E in physical space. As discussed in 0, this linking number is 
precisely I — the number of times a nontrivial fiber So twists around the loop's core, which 
has S = —In- Thus flux loops with an I = 1 planar twist form fundamental skyrmions; flux 
loops with nonplanar twist have / > 1, which is trivial for Nf > 3. 

A nicer parametrization of the skyrmion stems from the exact sequence 

^We consider here quarks whose trajectories remain well-separated from degeneracies at the origin; 
that is, for which r 2 >> m^ 2 , where rrii are the quark mass eigenvalues. Thus we rely on an 
asymptotic regime inside the flux tube, in which m~ l « m~ l . The flux tube background affects 
general quark wave functions more subtly. Since the mass eigenstates depend explicitly on 9, off- 
diagonal 9 ■ t terms appear in the quark equations of motion, mixing mass eigenstates at small 
radius. This causes corrections to maximal Aharonov-Bohm scattering, much as in Q. 
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n 3 (SO(N f )) - Tr 3 {SU{N f )) - n 3 (SU(N f )/SO(N f )) - n 2 (SO(N f )) = . 

This implies that fundamental skyrmions in the theory can be constructed from the funda- 
mental skyrmions 

gf(r, w) = exp (iF s (r) u { T in ) (5.2) 

in SU(Nf). Here r and u are the radius and unit direction vector in 3-space, and F s (r) 
approaches 2tt at r = and zero at r = oo. This determines an axisymmetric skyrmion 

E s = t - 2 sin 2 (F s /2) (1 - «*) 1 □ + 2i sin(F s /2) • 

(cos(F s /2) (u y T xn + u x r zn ) + sin(F s /2) u z (-u x r xU + u y r zD )) , (5.3) 

after a global spatial rotation fixing the z-axis as the axis of spatial axisymmetry. In the 
xy-plane, this gives 

Y, s (z = 0) = 11 + (cos Fg — 1) 1 □ + i sin F s (u y t xU + « x r zD ) . (5.4) 
Comparing with equation |2.11| for the flux tube allows us to identify the angular winding 



n&(#) Tf, n of a flux tube with that of a skyrmion in the rry-plane. Thus, if their radial 
boundary conditions coincided, we could deform the lower hemisphere of the skyrmion into 
a flux tube. However, because both F and F s vary from to 2it over this plane, the skyrmion 
cannot end in a single flux tube. The boundary conditions instead allow the skyrmion to join 
only to flux tube configurations where F(r) ranges from to 4tt — that is, configurations 
with two flux tubes, deformable to the trivial configuration. Thus skyrmions cannot be 
confined in this theory. 

However, objects which combine to form skyrmions can interact with the flux tubes. 
Such "half- skyrmions" could arise as external spinor sources in the underlying theory. They 
should be confined, as fundamentals cannot screen them. As mappings on G/H, they 
appear in our theory precisely as half-skyrmions, that is, as objects of the form [5]3] with 
F s (r) ranging from to 7r. Such objects are not defects in the conventional sense, since 
they have linearly divergent energy — just like an unscreened point source. Their boundary 
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conditions allow them to join to their opposite winding counterparts via single flux tubes, 
confining their linearly divergent energy to a length scale set by the tube length. We thus see 
that confinement of sources in an SO(N c ) gauge theory can induce a relic phenomenology, 
which persists in the low energy Skyrme model. 
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APPENDIX: 



In this appendix, we show that trivializations |2]3| of the loop h (a) produce minimal 



flux tubes only when h((3) = 11 . We do this by generalizing the analysis following equation 



for arbitrary h((3) — only to find that consistency and energy considerations restore the 
restriction h{f3) = 11 . 

Note that conjugation by h{(3) rotates 2Th in so(Nf): 

n (jk)(P) r y (j k ) -> u {jk) {(3) T y{jk) , where u {jk) ((3) = h( jp 0) n (pq) h {kq) (p). (Al) 

Since h{0) is real and orthogonal, uu k \{0) is just a real unit vector, like n^ k y b(0) then 
conjugates this rotated generator, giving Equation |2~5| with the substitution n(j k ) — > U(j k )(/3)- 
This fully conjugated generator must become — 2Th at j3 = tt, to produce the loop g{a,P = 
7r) = h 2 (a). This occurs only if uy k )(n) = \(yk) n (jk)7 f° r eacn (jk), with 



A 



— exp(il7i) for (jk) = (12) 

- exp (ilir/2) for j < 2, k > 3 (A2) 
-1 for 2 < j < k. 
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We note that these matching conditions depend only on a basis choice fixing b((3) in the 
(12) plane. Rotations of b(f3) within that plane, for example taking T b — > |r x ( 12 ), change the 
details of equation |2.5| , but not the coefficients of TyUjA and its orthogonal generators; and 
these coefficients alone fix the matching |X2 . 

These matching conditions cannot be achieved by orthogonal conjugation of an arbitrary 
generator T^. A first constraint stems from the reality of Xtjk)- This allows a nonvanishing 
n^jk) only when the relevant exponential gives ±1. For 1 = 1, this forces n^ik) = n^k) — for 
k > 2 — splitting both h(a) and h(/3) h(a) h _1 (P) into two commuting blocks. The 2 < j, k 
block is unaffected by b- conjugation, and changes g only by an overall right multiplication, 
which leaves S = gg T invariant. Thus, for I = 1, Th = |t"j/(12) again produces all nontrivial 
loops g. 

Equation [Al] places more subtle constraints on T h , related to consistency of the diagonal 
reduction u^k)^) = A^n^), for all {jk). This is solved, for independent rijk, only when 
^0»( 7r ) ^ s itself diagonal: h^ p )(Tt) = Aj5( JP ), with Xj real. This constrains the matching 
conditions attainable by nonzero to Xjk = XjXk- Note that our rijk are effectively 
independent, since the relation 1X21 holds for all global rotations of nu]A that leave b(f3) in 
the (12)— plane. Such rotations map a single generator in any of the three matching classes 

- (jk) coincident with (12);(j'/c) intersecting (12) in a single line; and (jk) disjoint from (12) 

- to the entire class, in arbitrary linear combinations. Thus, for I = 2, only two possibilities 
are consistent with the matching conditions [A2| . We can have nonvanishing n^^ in the cross 
planes (lk), (2k), with k > 3. However, this precludes nonvanishing nuk) outside the cross 
planes: X^) — A( 2 fc) = 1, for all k > 3, implies X(jk) = 1 for all j, k. Thus we must set 
ri(jk) = outside cross planes to obey |A2|. Similarly, nonzero n( 12 ), with A(i 2 ) = —1, implies 
Xfik) = — X(2k) for A; > 3 — forcing = in the cross planes to obey |X2|. The case 
/ = 2 thus produces no new candidates for T h : either nuk) vanishes outside the planes (jk) 
which intersect (12) in a single line, or n(ju) = e (i2) in the (12) plane (discarding a right 
multiplication g — > gh! as in the / = 1 case above). 

The only surviving consequence of h— conjugation, then, is an overall rotation U(jk)(n) = 

26 



±ri(jk)- For the positive sign, h((3) commutes with /i _1 (a:) at both endpoints in (3, so choosing 
[h((3),h^ 1 (a)] = for all (3 clearly minimizes gradient energy. For the negative sign, a 
similar reduction occurs for b{(3). For then b(/3) commutes with h(/3) h~ l (a) h~ l (j3) at both 
endpoints, and we can deform g to lie entirely in H. Thus only the case h(/3) = 1 induces 
a nontrivial defect £ G G/H of minimal energy. 
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FIGURES 

FIG. 1. a.) Flux tube solutions F(r) for A values 0.2, 0.236, 1.0, and 2.0, with mass stabilization 
only (uj = 0). The dotted line corresponds to the physical Ao = 0.236. Note that core size shrinks 
with increasing A. b.) The above flux tubes' energy density p. Note that tension grows with the 
coupling A. 

FIG. 2. Dependence on stabilizing mass (u = 0). Solid line shows variation of the total tension 
t = J d 2 r p as & function of A. Dotted line shows the gradient contribution To only. The physical 
Ao = 0.236 is marked. 

FIG. 3. Zero mode dependence on rotational direction. Solutions for A = Ao, with fixed 
frequency to = 1, and varying planar component n n . Note that core size shrinks and tension grows 
as the planar component n n increases from to 1. This implies that twisting dislocations are 
favored over slipping ones (see text). Allowing A to vanish has negligible effect on F and ptot- 

FIG. 4. The rotating mass-stabilized flux tube (A = 0.236). For n n = 0, shows the flux tube's 
deformation due to rotational zero modes. As for other values of A, the flux tube's core size shrinks 
and tension grows as the rotation frequency lo increases from 0.5 to 1.5. Rotation with a nonzero 
planar component n n displays the same trend, but begins with more compact and energetic flux 
tubes when uj = 0.5 (cf. Fig. 3). 

FIG. 5. Rotational deformation, for n n = 0.5. Dotted line shows Ao = 0.236; solid lines show 
the A = and A = 1 cases, a) Variation of the ground state tension r with rotational frequency 
uj. (Note that f, the tension conjugate to z of Section IV, is given by (eF 7r ) _1 r). b) and c) 
Variation of the planar and nonplanar moments of inertia A Q and A D /. (Note that A Q and A D / of 
section IV are given by 4(eF 7r )~ 3 times their counterparts plotted here.) 
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